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Рассматривается задача прогноза добычи нефти и оценки извлекаемых запасов с учетом дополнительной априорной информа- 
ции и экспертных оценок технологических параметров разработки месторождений. Приводятся результаты анализа точности 
оценок в зависимости от числа лет разработки. 


Введение 

В последние годы в нефтяной промышленности 
в связи с широким внедрением информационных 
систем, современных технологий исследования 
скважин, систем регистрации разнообразной ин- 
формации о состоянии объектов разработки, сква- 
жин, нефтяных пластов и месторождений возрос 
интерес к задачам построения математических мо- 
делей процессов нефтегазодобычи [1—6]. 

Проблемы построения математических моде- 
лей процессов нефтегазодобычи связаны с реше- 
нием задач идентификации, которые часто называ- 
ют обратными задачами. Задача идентификации 
заключается в построении оптимальных, в смысле 
заданных критериев качества, математических мо- 
делей технологических параметров разработки (до- 
бычи нефти, жидкости, воды, пластовых давлений, 
обводненности продукции и т. д.) с использовани- 
ем промысловых данных, результатов комплекс- 
ных исследований скважин и нефтяных пластов. 
Задачи идентификации технологических параме- 
тров подразделяют на два крупных направления, 
имеющие свои цели и задачи. 

Первое направление представляет задачи уровня 
проектирования разработки нефтяных месторож- 
дений, которые решаются крупными коллектива- 
ми в научных центрах нефтяных компаний, про- 
ектных институтах [6]. На стадии проектирования 
разработки месторождения создаются цифровые 
геолого-технологические модели нефтяных место- 
рождений, позволяющие прогнозировать показа- 
тели разработки на достаточно длительный период 
(20-30 лет), создавать технологические схемы и 
проекты разработки, определять стратегию разви- 
тия компании. 

Второе направление представляют задачи иден- 
тификации уровня мониторинга и оперативного 
управления разработкой нефтяных месторожде- 
ний. На стадии разработки месторождений акту- 
альным является моделирование технологических 
параметров для оперативного решения задач прог- 
нозирования добычи нефти, оценки эффективно- 
сти геолого-технических мероприятий, определе- 
ния оптимальных режимов работы скважин и т. п. 

Здесь более мобильными, легко адаптируемыми 
(настраиваемыми) на основе промысловых данных 
и результатов исследований скважин по сравнению 


с геолого-технологическими моделями уровня 
проектирования являются различные регрессион- 
ные статические и динамические модели техноло- 
гических параметров разработки. Наиболее широ- 
кое применение получили промыслово-технологи- 
ческие модели, основанные на характеристиках 
вытеснения, уравнениях фильтрации флюидов, ма- 
лопараметрических регрессионных моделях добы- 
чи нефти, жидкости и воды [7, 8]. 

Однако реальные условия нефтегазодобычи ха- 
рактеризуются неполнотой, неоднородностью, недо- 
статочностью исходных данных о фильтрационных 
параметрах и энергетическом состоянии нефтяных 
пластов, различными ошибками регистрации техно- 
логических и геологических параметров разработки, 
неадекватностью моделей и т. п. [1—5]. В данных 
условиях использование классических методов иден- 
тификации технологических параметров разработки 
часто связано с проблемами устойчивости и низкой 
точностью получаемых решений [5]. 

Перспективным направлением преодоления 
реально существующей сложности процессов неф- 
тегазодобычи, недостаточности данных, повыше- 
ния точности и устойчивости решений, является 
использование интегрированных систем иденти- 
фикации с учетом разного рода дополнительной 
априорной информации, накопленного опыта и 
знаний [9, 10]. 

Интегрированные системы идентификации 
предоставляют возможность объединять математи- 
ческие (физически содержательные) модели техно- 
логических параметров разработки месторождений 
и формализованные экспертные оценки лица, при- 
нимающего решения, в единую интегрированную 
систему моделей, отражающую целостные свой- 
ства реальных систем, что дает интегральный (си- 
нергетический) эффект при решении разнообраз- 
ных задач идентификации и управления [11]. 

Данная работа является развитием технологии 
интегрированных моделей при идентификации од- 
ного из важных технологических параметров разра- 
ботки месторождений - добычи нефти с использо- 
ванием малопараметрических регрессионных моде- 
лей отбора нефти и непараметрического метода 
аналогов [10] в условиях априорной неопределен- 
ности о статистических характеристиках дополни- 
тельных априорных данных и экспертных оценок. 
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Интегрированные модели 
и алгоритмы идентификации добычи нефти 

Рассмотрим интегрированную систему моделей 
добычи нефти на основе малопараметрической мо- 
дели, отражающей эволюцию текущего отбора 
нефти во времени, с учетом дополнительной ап- 
риорной информации об извлекаемых запасах, па- 
раметрах модели добычи нефти и прогнозных зна- 
чениях добычи, которую представим в удобном для 
исследований матричном виде: 

а ;= а »+^ 

_ т 

8 =з(Т,а)+г] = |/(т ,а)Ат +77 , 

( 1 ) 

аі =а +ѵ ; , / = 1 ...й, 

6„ =<2(т,а) + е. 


где і=\...п), (?„(а)=(/(/„а), і=\...п) - век- 

торы значений добычи нефти, полученные в про- 
цессе разработки месторождения и на основе моде- 
ли Д/, а) за соответствующие промежутки времени 
А /=/,—/,•_!, / = 1 ... /7 (год, месяц); 51 =Щ, ]=\...р), 
а,=(а р ;=1. ../77), 1=1. .У, 0=(б„(Г„+т 7 ), ]=\...г) - 
векторы экспертных оценок извлекаемых запасов 
нефти за время разработки Т, дополнительных ап- 
риорных данных о параметрах модели добычи 
а=(а ь а 2 ,...,а т ) и прогнозных значений добычи 
нефти, полученные из разных источников инфор- 
мации, например, из проектных технологических 
документов на разработку месторождений; 
%,г),ѵ,,1=\...сІ,е - векторы случайных величин с ну- 
левыми математическими ожиданиями, предста- 
вляющие погрешности измерения добычи нефти, 
ошибки, вызванные неадекватностью модели до- 
бычи нефти, ошибки дополнительных априорных 
сведений и экспертных оценок. 

Рассмотрим решение задачи параметрической 
идентификации модели (1) с использованием ком- 
бинированного критерия качества 

Ф(а,Р) = |е*-б„(а)|Г + |І5-5'(Г,а)|| + 

II II II иг, 

+^^ 2 Да/-а) 2 + Ь-/„(а)|| (2) 

і-і 


непараметрического метода аналогов [10] доста- 
вленного из частных квадратичных критериев ка- 
чества, где ||Х|р - норма X, 

Щ =Аш ё (К( $(Т ’ а ' п ] ~ $и ),]=1...р), 

Р\ 

а"„ - а я 

Щі = Лаі(^(^— , 7 = 1 . ..777), 

Рг 


Щ = с1іа<Ц/<'( 


+ т 1 ) 

А 


),./ = !••+)- 


диагональные матрицы весовых функций К(и) с 
управляющими параметрами /3,,/3 2 ,/3 3 , учитываю- 


щих значимость (вес) дополнительных априорных 
данных и экспертных оценок, а 0 =(а 1 0 ,а 2 0 ,...а^) - 
вектор начальных приближений параметров моде- 
ли добычи нефти. 

Функции (ядра) К(и) введены по аналогии с не- 
параметрическими оценками плотности вероятно- 
сти и регрессии, обладают следующими свойства- 
ми [12]: 

К(х-и)/ И) — >1, И — э оо; К(х— и)/ И) — э 0, А— >0. 

Задача параметрической идентификации (1) 
сводится к решению оптимизационной задачи по 
определению вектора параметров модели добычи 
нефти 

а\р) = аг§тіпФ(а,/1) (3) 

а 

и оптимизационной задачи по определению векто- 
ра управляющих параметров 

Р* =аг § тт(|е*-б„(а*,/3)|р), (4) 

где ащ тіп/(х) обозначает точку минимума х функ- 
ции Дх). 1 

Приведем решение оптимизационной задачи 
(3) с использованием метода Ньютона [13]: 

а' =а‘~' + А ; Аа‘ _1 , і =1,2,3,..., (5) 

где оценки приращений вектора параметров ДсГ 1 
на каждой итерации определяются путем решения 
системы линейных уравнений 

А'-'Аа 1 -' = В'-', (6) 


где 


А 1 -' = (О т н О н + ГУ ГД + X + ГУ IV, о _ у- 

>1 

В 1 -' = (В[е н +е,ІѴД + ^ ѵ Аа, + ё г - ГД У~\ 


О = 


д/(у,а) 

V да і 


І = 1.../7, 7 = 1. „777 


Д = 


д/(1 п +т п а) 

да, 


І = 1...Г,7 = 1 .../77 


О = 


матрицы 

д$(Т,а) . 
7 = 1...777 

да, 


частных производных; 

Л 

- вектор частных произ- 

V 


водных от функции извлекаемых запасов $(Т,а); 
е„=( ё=(5-з(Т,а 0 )), Аа,=(а,—а°), 

Ѵ~(0 іГ 0( г, а")) - векторы невязок. 

Для получения системы линейных уравнений 
(6) необходимо: 

1. Разложить функции Д?, а), $(Т,а) в ряд Тейлора 
и ограничиться членами первого порядка мало- 
сти относительно приращения вектора параме- 
тров Да 0 : 
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/(/,а) = Д/,а°) + ]Г(^^) 0 Да°, 

у '=1 

5(Г,«) = 5(Г,«°) + ^( &(Г ’ а) )°Аа°. (7) 

7^ За. 

2. Используя разложение (7) перейти от модели 
(1) к линейной интегрированной системе моде- 
лей добычи нефти относительно приращения 
вектора параметров Да 0 : 

е я =77 й °Да 0 + %, 
еі = Г>“Да° +т], 

■ -о (8) 

Да/ =Да +Ѵ[,1 = 1...Д, 

ё° = О 0 Да° + е. 


3. Для модели (8) по аналоги с (2) составить ком- 
бинированный функционал качества 

І|2 


Ф(Да и ) = ||е н °-2)„ 0 Да 011 


е, - Д Да 


+^Г 2 ,(Да°-Да°) 2 + е г °-Д°Да' 


2 

1Г, 


(9) 


4. Вычислить частные производные от функцио- 
нала (9) по параметрам Да 0 и приравнять их к 
нулю 


дФ(Ла°) 

(ЭДа° 


-(0° Н ) Т УѴ\(е° н -О°Аа 0 ) - 


-(оуиг (е, - о" Да)) - ^ Г 2 , (Да “ - Д а 0 ) - 
/=1 

-(Д 0 ) г Г 3 (е° г -Д 0 Да 0 ) = 0. (10) 


После проведения несложных преобразований 
из (10) следует система линейных уравнений (6). 

Следует отметить, что для решения оптимиза- 
ционной задачи (4) по определению вектора упра- 
вляющих параметров целесообразно использовать 
методы оптимизации без вычисления производных 

[13]. 


Точность оценок прогноза добычи нефти 
и извлекаемых запасов 

Рассмотрим вопросы точности алгоритмов 
прогнозирования добычи нефти и оценки извле- 
каемых запасов на примере интегрированной сто- 
хастической системы моделей (1), которую пред- 
ставим в удобном для статистического моделирова- 
ния матричном виде: 

б* = б„ (а )(1 + с^), 

8 = з(т,а)(1 + с 2 гі), 

- , , (П) 

а/ = а(1 + с,)ѵ,,/ = 1...а, 

Я„ =б(г,а)(1 + с 4 )е, 

где 0> и ,8 ,а„!=\ ...сІЯ„ - векторы имитируемых значе- 
ний добьии нефти, дополнительных априорных дан- 


ных и экспертных оценок извлекаемых запасов, пара- 
метров модели и прогнозных значений добычи, полу- 
ченные с использованием функции регрессии [7, 10], 

/(Г, а) = а 1 ехр(-а 2 /)-/ С(з , (12) 

широко используемой при моделировании промы- 
словых данных в качестве представления зависи- 
мости годовой добычи нефти от времени разработ- 
ки; — векторы случайных величин, 

полученные с использованием датчика псевдослу- 
чайных чисел, распределенных по нормальному за- 
кону с нулевыми средними значениями и единич- 
ной дисперсией; с р /-1...4 константы, предста- 
вляющие относительный уровень ошибок. 

Для определения точности оценок прогноза до- 
бычи нефти 

&*(/„+ г) = /(/„+ г, а*(/Г)) (13) 


и извлекаемых запасов 


Д* =.[/(/, а*(Д*))Л 

О 


(14) 


использованы их относительные ошибки 


Д/„+т,а) 


5„(Д*) = аЬ 8 ( 


з(Т,а) 


(16) 


в зависимости от длительности разработки 1=п (го- 
ды) при значениях параметров добычи нефти 
а!=50000, а 2 =2,1, а 3 =0,26. Здесь а*(/Г) - оценки па- 
раметров модели добычи нефти (12), полученные по 
алгоритму (6). Для определения оптимальных значе- 
ний управляющих параметров /Г (4) использован 
метод деформированного многогранника [13]. 

В табл. 1 приведены результаты расчета относи- 
тельных ошибок прогноза добычи нефти ( 15) ( т=3) и 
извлекаемых запасов (16) в зависимости от вида ин- 
тегрированной модели добычи нефти (11) и уровня 
ошибок дополнительных априорных сведений. Для 
указания вида модели добычи нефти введены сле- 
дующие обозначения: { б„*} - модель добьии нефти 
(априорная информация не учитывается); 
{0' н ,8,а}, ~ интегрированные системы 

моделей с учетом экспертных оценок извлекаемых 
запасов, дополнительных данных о параметрах мо- 
делей и прогнозных значениях добычи нефти, для 
которых введены условные обозначения М1-М4. 

Для оценок параметров а модели добычи неф- 
ти М1 и начальных приближений параметров а 0 в 
алгоритме (5) использован метод Гаусса-Ньютона 
[14] 

а' = а'~‘ + /г,.Да ,_1 , /=1,2,3,..., (17) 


где на каждом шаге / решается система линейных 
уравнений вида 

(0:0 н Г-Аа ! -'=(0 т н е н У-\ 

которая следует из (6) при /3,=0,/=1...3. 
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Уровни ошибок дополнительных априорных 
данных и экспертных оценок имитировались сле- 
дующие: добыча нефти 5 % (^=0,05), извлекаемые 
запасы 2, 5 и 10 % (с 2 =0,02; 0,05; 0,1), дополнитель- 
ные априорные данные параметров модели добычи 
-2 % (с 3 =0,02), а ошибка прогноза добычи нефти 
равна 5 %, что соответствует приемлемым на прак- 
тике уровням точности оценок. Длительность раз- 
работки составляла три года (п= 3). Количество до- 
полнительных данных и экспертных оценок в (11) 
для каждого типа данных равнялось трем 
(р=ё=г= 3). Вес дополнительных априорных дан- 
ных и экспертных оценок при формировании ма- 
триц Щ, ]Ѵ 2 1 , Щ (6) рассчитывался с использовани- 
ем функции К(х)=ехр(-х 2 ). 

В табл. 2, 3 приведены значения относительных 
ошибок прогноза добычи нефти ( 1 3) (г=3) и извле- 
каемых запасов (14) в зависимости от числа лет 
разработки и вида интегрированной модели добы- 
чи нефти. Уровень ошибок имитируемых значений 
дополнительных данных об извлекаемых запасах и 
прогнозных значениях добычи составлял 5 %. Уро- 
вень ошибок имитируемых значений дополнитель- 
ных априорных сведений о параметрах модели до- 
бычи нефти составлял 2 %. 


Таблица 1. Относительные ошибки оценок прогноза добычи 
нефти и извлекаемых запасов, % 


Оценки 

Вид интегрированной модели добычи 
нефти 

1 

2 

3 

4 

М1 

М2 

М3 

М4 


Уровень ошибок экспертных оценок 
извлекаемых запасов 

2 

5 

10 

2 

5 

10 

2 

5 

10 

Прогноз добычи 

53 

2,6 

3,2 

5,1 

2,2 

2,8 

4,5 

1,9 

2,3 

3,1 

Извлекаемые запасы 

49 

2,3 

3,1 

4.6 

1,8 

2,6 

3,5 

1,4 

2,2 

2,8 


Таблица 2. Относительная ошибка оценок прогноза добычи 
нефти, % 


Модели добычи 
нефти 

Длительность разработки (номер года) | 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

М1 

- 

- 

ш 

ш 

ш 

16,1 

11,3 

8,7 

6,9 

5,3 

М2 

ш 

ш 



гаи 

9,7 

7,1 

5,6 

4,8 

4,2 

М3 

ш 

ш 




7,4 

5,2 

4,7 

4,1 

3,6 

М4 

ш 

ш 



ш 

3,9 

3,1 

2,6 

2,3 

2,1 


Таблица 3. Относительная ошибка оценок извлекаемых запа- 
сов, % 


Модели добычи 
нефти 

Длительность разработки (номер года) | 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

М1 

- 

- 

49 

37 

21 

14,6 

9,3 

7,2 

5,7 

4,9 

М2 

46 

34 

20 

15,6 

11,6 

9,5 

7,2 

5,9 

4,6 

3,8 

М3 

40 

20 

17 

12,3 

10,2 

8,1 

5,9 

5,2 

4,6 

3,3 

М4 

25 

13 

8,8 

6,5 

5,1 

3,7 

3,2 

2,8 

2,3 

1,9 


Из табл. 1-3 видно, что учет дополнительных 
априорных данных и экспертных оценок приводит 
к существенному увеличению точности оценок 


прогноза добычи нефти и извлекаемых запасов на 
ранних стадиях разработки месторождения 
(3-6 лет) по сравнению с оценками (17) на основе 
модели добычи нефти М1, где априорная инфор- 
мация не учитывается. Например, относительная 
ошибка извлекаемых запасов, полученная к третье- 
му году разработки с использованием метода Гаус- 
са-Ньютона (17), составляет 49 % (табл. 3), а с ис- 
пользованием метода интегрированных моделей с 
учетом априорной информации (5) порядка 9 %. 

Приведем результаты прогноза добычи нефти с 
использованием фактических промысловых данных 
годовой добычи нефти одного из месторождений 
Тюменской области за 15 лет его разработки с 1990 г., 
изображенных на рисунке (линия 1). Линия 2 пред- 
ставляет прогнозные значения добычи нефти (13) 
по данным трех лет разработки п = 3 и т=1,2,...,12, 
полученные с использованием модели (12) и алго- 
ритма (17). Линия 3 отражает прогноз добычи неф- 
ти, с использованием дополнительной априорной 
информации об извлекаемых запасах 6,3 млн т, рас- 
считанных до начала разработки месторождения. 



Номер года разработки 

Рисунок. Фактические значения добычи нефти (линия 1), 
прогнозные значения добычи с использованием 
модели МТ и М2 (линии 2, 3) 

В табл. 4 отражены значения относительных 
ошибок прогноза добычи нефти, начиная с четвер- 
того года разработки месторождения, полученные 
с использованием модели М1, алгоритма (17) и мо- 
дели добычи нефти с учетом априорной информа- 
ции об извлекаемых запасах М2 и алгоритма (5). 


Таблица 4. Относительная ошибка оценок прогноза добычи 
нефти, % 


Модели добы- 
чи нефти 

Длительность разработки (номер года) 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

М1 

І4 

11 

9 

8,5 

7,4 

4 

Т 

2,9 

2,5 

1,8 

1,7 

1,6 

М2 

~6~ 

4,4 

3,8 

4,1 

3,6 

2,5 

18 

2,1 

1,9 

1,1 

0,9 

0,8 


Из рисунка и табл. 4 видно, что использование 
дополнительных априорных данных об извлекае- 
мых запасах позволяет существенно повысить точ- 
ность прогноза добычи нефти на ранней (первой) 
стадии разработки месторождения. 
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Технические науки 


Выводы 

Предложенные интегрированные модели добы- 
чи нефти и алгоритмы идентификации позволяют: 

1) учитывать дополнительные априорные данные 
и экспертные оценки технологических параме- 
тров разработки месторождений: извлекаемые 
запасы, прогнозные значения добычи нефти, 
параметры модели добычи нефти и т. и.; 

2) получать оценки прогноза добычи нефти и из- 
влекаемых запасов в условиях априорной нео- 
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